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SPECIFICATION 





This method of reducing multipath avoids the cross term noise produced by other 
methods [Refs. 1 through 10]. These other methods correlate (cross correlate or auto 
correlate) two versions of the signal each containing the direct path and multipath 
reflections. A correlation peak between the direct path and a reflection gives a 
measurement of the relative delay between the direct path and the reflected path and the 
amplitude and phase of the reflection. As the relative delay between the two signals is 
moved from this peak, other correlation peaks are produced when there is correlation 
between the reflections, i.e., when the relative delay equals the difference of the delays 
for two objects. This produces noise which interferes with the measurement since there 
are many possible combinations/ When one delays one of the signals and subtracts it 
from the other with the correct amplitude and phase, he tends to cancel the multipath. 
The other reflections, however, are dragged along with it and produce noise. These 
noises are avoided by the invention described in this disclosure. 



The Appendix contains listings of computer programs written in the APL language and 
a sample computer run. The programs embody the basic features of this invention. The 
computer run is a demonstration of its performance on simulated data. 

The program NORM generates simulated data. The program DELRES calculates the 
delay residuals. The program NLS minimizes the residuals defined in DELRES by 
nonlinear least squares. 

The program NLS is first edited so as to output the residuals as RR. The program 
NORM is used to generate 100000 mean zero variance 1 random numbers in X5. This 
represents a 50000 Hz broadband signal sampled at 100000 samples per second with 
r.m.s. (root mean square) of 1. Two multipaths are simulated both of 0. 1 magnitude. One 
is delayed 1 sample; the other is delayed 3 samples. The result is put into XX5. The 
signal to multipath noise ratio is calculated to be 17 dB. The NLS program is applied to 
XX5 and the two multipath noises are canceled to a level of 51 dB signal to multipath 
noise ratio. 



Statement Regarding Federally Sponsored R&D RECEIVED 

No Federal funds were used to conceive or develop this invention. 
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Background 

The direct path of a radio signal from transmitter to receiver is frequently interfered 
with by reflections of the signal from stationary and moving objects. This is called 
multipath noise. This invention utilizes a new adaptive filter technique to reduce 
multipath noise. This method of reducing multipath avoids the cross term noise produced 
by other methods [Refs. 1 through 10]. These other methods correlate (cross correlate or 
auto correlate) two versions of the signal each containing the direct path and multipath 
reflections. A correlation peak between the direct path and a reflection gives a 
measurement of the relative delay between the direct path and the reflected path and the 
amplitude and phase of the reflection. As the relative delay between the two signals is 
moved from this peak, other correlation peaks are produced when there is correlation 
between the reflections. This produces noise which interferes with the measurement. 
When one delays one of the signals and subtracts it from the other with the correct 
amplitude and phase, he tends to cancel the multipath. The other reflections, however, 
are dragged along with it and produce noise. These noises are avoided by the invention 
described in this disclosure. 

Brief Summary 

A nonlinear least squares method measures the delay, Doppler shift and amplitude of 
the multipath due to each object and subtracts a very accurate reconstruction of each 
multipath signal from the noisy signal. If an object is a target, its range, range rate and 
magnitude is got from the measured multipath delay, Doppler shift and amplitude. 
Position and velocity of the target can be obtained by geometric triangulation with 
multiple transmitters. Target angle can be measured by the relative phase of the 
corresponding filter coefficients across multiple antennas. The system can be used on a 
surveillance aircraft to cancel ground reflections and measure targets. 



Brief Description of Drawing 

The Figure is a block diagram illustrating the adaptive filter of the invention to reduce 
multipath in a radio transmitted signal wherein the multipath is caused by reflections 
from fixed and moving objects, interferers, or targets. 

The method works best on white signals; therefore, prewhitening and post unwhitening 
filters are employed. 

The differential delay and Doppler shift for each target is measured by the filter weight 
x mn which corresponds to the differential delay nDand Doppler shift mf caused by that 
target. The magnitude of the filter weight is a measure of the strength of the target. 
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The quantity labeled "Residual" in the Figure is the received signal minus the result of 
application of the filter weights w l9 w T9 w 3 . 9 x U9 x l29 x l y-' 9 x 2l9 x 22 . 9 x 2 ^'-x 3l9 Xy 29 Xy i '-'' . A 
nonlinear least squares method is used to pick the filter weights which minimize the 
mean squares of the "Residual". The result is a clean version of the original signal which 
is the received signal minus the multipath reflections. 

The filter weight x^ corresponding to a moving target is a measure of the differential 
range nDc and range rate mfc where c is the speed of light. The system, with a 
common antenna, can be duplicated to receive two signals each from each of two widely 
separated transmitters. Geometric triangulation can be used to measure the two 
dimensional position and velocity of the target. Triplication can be used for three 
dimensional position and velocity. 

Another use of the system is to detect the angle of moving targets. If multiple antennas 
are provided, each connected with a system like that shown in the Figure, the angle of a 
moving target causing a multipath reflection to the receiving antennas can be measured 
by using the relative phase of the corresponding delay and Doppler complex coefficients 
across the several antennas . 

The receiving antennas, receivers and processing system can be placed in a surveillance 
aircraft the position and velocity of which is obtained by an accurate navigation system 
such as GPS. Objects on the ground, interferers and targets causing multipath reflections 
can be processed by the system. Target position and velocity can be obtained by adding 
the position and velocity of the surveillance aircraft to the measured position and velocity 
of the target. 



Detailed Description 

This method of reducing multipath avoids the cross term noise produced by other 
methods [Refs. 1 through 10]. These other methods correlate (cross correlate or auto 
correlate) two versions of the signal each containing the direct path and multipath 
reflections A correlation peak between the direct path and a reflection gives a 
measurement of the relative delay between the direct path and the reflected path and the 
amplitude and phase of the reflection. As the relative delay between the two signals is 
moved from this peak, other correlation peaks are produced when there is correlation 
between the reflections This produces noise which interferes with the measurement. 
When one delays one of the signals and subtracts it from the other with the correct 
amplitude and phase, he tends to cancel the multipath. The other reflections, however, 
are dragged along with it and produce noise. These correlation noises are avoided by the 
invention described in this disclosure. 
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The Figure is a block diagram illustrating the adaptive filter of the invention to reduce 
multipath in a radio transmitted signal wherein the multipath is caused by reflections 
from fixed and moving objects, interferes, or targets. 

The method works best on white signals; therefore, prewhitening and post unwhitening 
filters are employed. 

The differential delay and Doppler shift for each target is measured by the filter weight 
x mn which corresponds to the differential delay nD and Doppler shift mf caused by that 
target. The magnitude of the filter weight is a measure of the strength of the target. 

The quantity labeled "Residual" in the Figure is the received signal minus the result of 
application of the filter weights W\?W2* w i* x \\* x \2* x n , '^ x 2\> x 72* x 2s"' x i\> x i2> x K"* * A 
nonlinear least squares method is used to pick the filter weights which minimize the 
mean squares of the "Residual". The result is a clean version of the original signal which 
is the received signal minus the multipath reflections. 

Refer to the Figure. The signal is modulated and amplified by the transmitter, Tx [Ref. 
17, ROHDE, p. 14]. After transmission from the antenna most of the power takes the 
direct path to the receive antenna. Some power, however, is reflected from stationary 
and moving objects such as buildings and aircraft. This tends to add or subtract from the 
direct path causing what is known as multipath interference noise[Ref 18, SKOLNIK, p. 
18.46]. This invention removes the multipath noise from the received signal. 

The signal at the receive antenna is amplified and demodulated by the receiver, Rx 
[Ref. 17, ROHDE, p. 361 ][Ref. 19, VAN VALKENBURG, p. 23-6}. It is then passed 
through a presampling filter, sampled and converted from analog to digital [Ref. 19, 
VAN VALKENBURG, p: 28-7]. These real digital samples have a Fourier transform the 
real part of which is symmetric about zero frequency and the imaginary part of which is 
antisymmetric about zero frequency. The real samples are first frequency shifted by 
multiplying by e~ J2nFt where t is time, F = and y is the sampling frequency [Ref. 

14, PAPOULIS, p. 15]. These complex samples are then passed through a low pass 
digital filter with a maximum frequency of F [Ref. 16, RABINER, p. 100][Ref. 13, 
OPPENHEIM, p 261]. The output of the low pass filter is resampled at a frequency of 
y D , i.e., every other sample, to obtain the complex samples of the received signal. 
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The received signal, after demodulation, analog to digital conversion and conversion to 
complex samples, can be modeled as t 

n=\ m=- hi 

where 

s t - original signal 

/ = sample index 

D = sampling time interval 

n = time delay index 

f = Doppler shift frequency interval 

m= Doppler shift index 

- Multipath amplitude coefficient 
j- unit imaginary number 

e = base of natural logarithms 
n= pi 

E = e n * 
t = iD 

7] t = white noise 
From the Figure, the residual is 

L N M 

where the w t are whitening filter coefficients and the x m are the delay-Doppler filter 
coefficients. This equation corresponds to that of Box [Ref. 1 1, BOX, p. 498]. One can 
minimize the root mean squared (r.m.s.) residual 

x)ver 

W \i W 2i"*<> W L > X -M\* X -M2>'"-> X -MN > X -M+U> X -M+\2 >' " -> X -M+\N " > X A4N 

by nonliear least squares [Ref 11, BOX, p. 504, p. 505] [Ref 12, MARQUARDT, p. 
431]. 
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To show that this removes the Multipath and restores the original signal, define the 
whitened signal as 

L 

Pi = s i ~H w i s i-i 

l=\ 

and the modified noise as 

L 

ft = ^- S w / l »-/ 

i=\ 

Then, the residual is 

1 =/»,+! Z °-e- J3m¥ "P,-.E m + <;,-t iv^ m 

n=\ m=~M n=l m=-M 

The root mean square residual is minimized when 

Jt » ae- j2mfhD 

nvi TTOl 

and 

r, * a + £ 

Note that only delayed samples of the whitened signal p i , i.e., p x _ n , /? * 0 , are 
subtracted from p. , so that it is not possible to cancel the whitened signal p t since it is 
uncorrelated with delayed versions of itself. 

Adding the above two equations which describe the action of the whitening filter, one 
obtains 

L 

1 * ( s r + - X w l( s r-l + *&-/)■ 
/=i 

Solving this equation, one gets 

which, from the Figure, gives the output of the unwhitening filter. This is the desired 
result which restores the original signal without multipath. The imply target 
differential range, range rate and magnitude. 



j 



6 




A simple test of this concept was made with the computer program shown in the 
Appendix. The Marquardt method for nonlinear least squares [Ref. 1 1-, BOX, p: 504, p. 
505] [Ref. 12, MARQUARDT, p. 431] was programmed in the APL language in program 
NLS. The equation for the residual was programmed in DELRES with 
r = 0 when m*0 

and 
w, =0 
If 

w l * 0 

this corresponds to the equation of Box [Ref 11, BOX, p. 498] and is a general recursive 
digital filter [Ref. 15, PAPOULIS, p. 45]. 

The original signal s i was generated with a Gaussian random number generator in 
program NORM with root mean square (r.m.s ) of 1 and / = 100000 . This corresponds to 
one second of a 50000 Hz bandwidth signal represented by 100000 real samples. The 
data is white, zero mean, Gaussian with r.m.s. = 1. 

The received signal q t was generated with 
a o\ = a oy - 0.1 , otherwise = 0 
and 

This corresponds to two multipaths both of magnitude 0.1, one is delayed 1 sample, the 
other is delayed 3 samples. The white signal itself provides an effective noise which 
perturbs the estimates. 

The process converged in three iterations, the final r.m.s. residual = 1.002370568. The 
three coefficients x 01 9 x m , jc 03 are shown to be approximately = 0. 1, 0.0, 0. 1 matching 
a m ,a 12 ,a 03 with r.m.s. errors of about 0.003. The mean residual was about 0.0087. The 
original signal to multipath noise ratio was 17 dB. After processing, the two multipaths 
were canceled to a signal to multipath noise ratio of 51 dB. 

The filter weight x^ corresponding to a moving target is a measure of the differential 
range nDc and range rate mfc where c is the speed of light. The system, with a 
common antenna, can be duplicated to receive two signals each from each of two widely 
separated transmitters. Geometric triangulation [Ref. 18, SKOLNDC, p 25.5, p 25.13] can 
be used to measure the two dimensional position and velocity of the target. Triplication 
can be used for three dimensional position and velocity. 

Another use of the system is to detect the angle of moving targets. If multiple antennas 
are provided, each connected with a system like that shown in the Figure, the angle of a 
moving target causing a multipath reflection to the receiving antennas can be measured 
by using the relative phases of the corresponding delay and Doppler complex coefficients 
across the several antennas [Ref. 18, SKOLNDC, p.3.34, p.3.35]. 
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The receiving antennas, receivers and processing system can be placed in a surveillance 
aircraft the position and velocity of which is obtained by an accurate navigation system 
such as GPS. Objects on the ground, interferers and targets causing multipath reflections 
can be processed by the system. Target position and velocity can be obtained by adding 
the position and velocity of the surveillance aircraft to the measured position and velocity 
of the target. 
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APPENDIX 
Computer Programs and Sample Run 



VN0BMLD1V 

[0] X+NORM NxZxTxR 

[1] p ;; 

[2] A 961204.1712 

[3] a GENERATE N NORMAL RANDOM VARIABLES, MEAN 0, SIGMA 1 

[4] a URL IS SEED 
[5] Z+2,\N*2 

[6] Z+Zp( 410000 )x? ( x/Z)p 10000 

[7] T+(o2 )x£[l; ] 

[8] -2x©£[2; ] )*0. 5 

[9] jsr^+,(i io,xj?)x2 lo.or 

V 1996-12-10 11.00.01 (G#!F-4) 

1DELRESIU11 
[0] DELRES X;N;I 

CI] A ; ; 

[2] P 000323.2241 

[3] q DMiUr RESIDUALS 

[5] i?^-oxr 

[6] Rl\Nl+Yl\Nl 

C 7 ] i>0 

[8] NEXTIi J«-J+l 

[9] i?[^+i]^-rcj\r+i]-^[^+r- ijn + .x* 

[10] +NEXTI*\I<( pY)-N 
[11] 

7 2000-03-25 22. 51.47 (GMT-^) 

[0] XE+XD NLS RESiX;bXxR;P;J;&RiA;B;DiEMiE;IiEEiAA 

[I] p iRES; 

[2] A 800225.0914 800820 000321.1532 

[3] fl NONLINEAR LEAST SQUARES 

[4] q CfMi?. RES IS RESIDUAL FUNCTION NAME 

[5] X+XDHil 

[6] AAVjri>[2;D 

[ 7 ] 1 ' 

[ 8 ] Jh-0 

[9] NBXTIiI+I+1 

[10] R+*RES 9 ' JST' 

[II] P (i? / BH*i?£S, ' A*)' 

[12] P EE+EE, (4/ + /(J?*2),[0.5]l)*0.5 

[13] (*/ + /(£* 2), [0.5]1)*0. 5 

[14] a AA^LyTJC RESIDUAL PARTI ALS 

[15] A P+*(RES 9 ? P f ), ' AT' 

[16] A+4JWMJ? 

[17] a CALCULATE R 9 A AND B 

[18] a (R A B)+±RES 9 ' (X LX ) 1 

[19] a +CALCRAB 

[20] A NUMERIC RESIDUAL PARTI ALS 

[21] F«- ' 1 

[22] ef+0 

[23] NEXTJzJ+J+l 

[24] hR+(*.RES 9 ' X+LXUlxJ=\f>X> )-R 
[25] P+P t Ai?*AJr[cn 
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t 



[26] +NEXTJ*\J<pX 

[27] p^((pjn,(pP)*pjn P p 

[28] ANPAR:A+(ISiP)+.xP 
[29] B*-($P)+.x/? 

[30] CALCRAB: 

[31] D+*(((0 l + pil)pi4)[;l] + l£-10)*0.5 

[32] AA+(A*Do .xD)+lE-10x(pA)pl 9 0*A 

[33] JKVJ!f-0x(DxB)S^ 

[34] ^^2 T Jx l J<3 

[35] R+slRES, » A-' 

[36] R (i? 4 B)+$lRES 9 *{X LX)* 

[37] £AK*/ + /1?,[0.5]1 

[38] £<-(*/ + /(J?*2),[0.5]l )*0.5 

[39] *^(JSr,£Af),[0.5]£x( ( (0 l + pAA)p(D° .xD)x$AA)L ;1],1)*0.5 

[40] P 
[41] £ 

[42] R » 1 

[43] P JST£f 

[44] p +NEXTI*\I<20 

[45] P 

7 2000-03-29 23 .00 . 55 (<W!r-4) 

[45] P RR+R 
[45] J?2?^-i? 

[46] 7 



pXS+NORM 1000 00 

100000 

□«-$+(*/ + /( (*5 )* 2 ),[. 5 ]1)*. 5 
1.002357739 

p^5^5[3+i99997] + .lxjf5[( x99997)o.+-l + 1 i+]+.xl 0 10 

99997 

□*jr^(*/ + /((JTJSr5-3+A'5)*2) > [ .5]1)*.5 
0.1418755036 

20xl0©5xJ7 
16.98230664 

(2 3p0 0 0 .001 .001 .OQDNLS » XX5 DELRES' 
1 .011856929 
1.002412348 
1.00237057 
1.002370568 

0.09896411424 0.001084301862 0.09780 701046 0.008745036907 
0.003147221067 0.003163105505 0.0031469113 59 1.002370568 

[WIM */ + /( (J?J?- 6 +JT5 )*2),[.5]1)*.5 
0.002841730873 

20x10®5tJV 
50.94879605 



